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ABSTRACT: Ticks have an exceptional ability to parasitize diverse animal and human lives and cause direct losses by transmitting 
a wide variety of pathogens to their host. For the management of the tick a vaccine using mid-gut antigen, Bm86 was developed 
which provided a variable degree of efficacy against different strains of the cattle tick, Rhipicephalus microplus. Contrary toa 
single antigen-based vaccine, it is expedient to use multi-epitopes vaccine targeting multiple functions of ticks. To develop a 
universal tick vaccine against multi-species, multiepitopic construct using immunodominant conserved tick proteins with B 
and T-cell epitopes were designed employing various immuno-informatics techniques. The 3D structure of the final construct 
was predicted, refined, and validated by different computational techniques. Molecular docking of the construct was carried out 
with TLR-2 receptor to predict efficient cell mediate the immune response. The findings manifest that the multi-epitopic 
construct might be active against multi-tick species infestations, a problem faced by many tropical and subtropical countries. 
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INTRODUCTION 

In the wake of modernization, endless human-driven 
changes in the microcosm and macrocosm have catalyzed 
the emergence, re-emergence, and spread of a myriad of 
new zoonotic diseases. As an analogy, nearly 30 
arthropod-borne human diseases are well recognized and 
many of them have emerged over the past years (WHO 
2020). Quite often, most such pathogens propagate 
unnoticed in enzootic tick-vertebrate cycles in wild natural 
conditions but may pose significant morbidity and mortality 
when spilling over to domestic animals or humans (Jahfari 
and Sprong 2016). Undoubtedly, ticks and tick-borne 
diseases pose a major “One Health” issue for humans 
and animals, worldwide, as they can disseminate an 
extensive array of pathogens causing anaplasmosis, 
babesiosis, Lyme disease borreliosis, Rocky Mountain 
spotted fever, Crimean-Congo hemorrhagic fever, 
Kyasanur forest disease, Tickborne encephalitis (TBE), 
etc. and incite direct damage to its host (Kumar et al. 


2020, Zhao et al. 2020). However, the concept of one 
health has not been suitably implemented due to 
insufficiency or the absence of transdisciplinary 
coordination in many developing countries. The problem 
is handled independently by animal and human health 
professionals and thus the magnitude of the problem has 
been further multiplied and reached a very alarming stage 
were at least 70% of emerging and re-emerging diseases 
are either zoonotic (spread from animals to humans) or 
vector-borne (carried from infected animals to others 
through insects). 

However, single recombinant antigens are seldom 
sufficient against complex parasites like ticks. The Bm86 
antigen, from the gut epithelium of R. microplus showed 
an appreciable but variable degree of efficacy (0 to 91%) 
in various geographical areas (de la Fuente et al. 2007). 
Practically, almost all single antigen vaccine using native 
antigen and their recombinant form is inefficacious to 
combat the infection (Willadsen 2008). Hence, contrary 
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to a whole (single) antigen-based vaccine, it is expedient 
to use multi-epitope vaccines that have considerable 
advantages such as safety, specificity, and cost- 
effectiveness. 

To address this affliction, the implication of protein 
derived subunit vaccine that includes short but remarkably 
immunogenic peptides seems promising as ticks are too 
complex to be used as a whole vaccine. Although 
pioneering technologies and a huge database are available, 
it is yet to be fully explored in tick vaccine development. 
A comprehensive understanding of tick immunity might 
give a leg up to developing highly successful vaccines. 
For instance, TGE (tick salivary gland extract) or saliva 
is reported to have immune-modiulatory, anti-hemostatic, 
and vasodilatory bioactive compounds possessing potent 
pharmacological activity which are released to temper 
and exempt host immune response for their co-survival 
with the host (Bowman ef al. 1997). Surprisingly, these 
chemicals subdue the pro-inflammatory response of 
macrophages by suppressing cytokine release and 
dramatically alter the proliferation of B and T lymphocytes 
(Ferreira and Silva 1998). Moreover, multiple components 
from the saliva of different ticks have been found to 
modulate host natural killer (NK) cells and complement 
system activation (Oliveira et al. 2010). Moreover, it has 
been recorded that dismantling the regular migration, 
maturation, and differentiation of dendritic cells (DCs) in 
vitro/in vivo by inhibiting IL-1, TNF-s< and by propagating 
Th-2 polarised state is characterized by pronounced IL-4 
and IL-10 level and diminished IFN-y and IL-12 level 
(Hannier et al. 2003). The DCs are proficient antigen- 
presenting cells and chief enablers of T-cell activation 
thus, their modulation provided an appropriate opportunity 
to ticks for evading the host’s immune responses that settle 
up with the host for the unimpeded meal and pathogen 
transmission (Cavassani et al. 2005). On this account, 
any halt in the precedent accord between the host and 
tick could result in a feasible tick control strategy. Here 
we propose to design multiepitope peptides against multi- 
tick species using CD8+ T cell epitope in combination 
with B and HTL epitopes from different antigenic proteins 
in an attempt to subvert this immunomodulatory state to 
the pro-inflammatory state. The multiepitopic peptides 


predictably have the potential to elicit strong cellular and 
humoral responses in hosts to combat multi-tick species 
infestation. Very similar initiatives pursuing a strategic 
approach to discover the epitope ensemble vaccine against 
SARS-CoV-2019 (Bibi et al. 2021), Ebola virus (Ullah et 
al. 2020), Acinetobacter baumannii (Solanki and Tiwari 
2018), Leishmania, Onchocerciasis and Plasmodium 
falciparum have been reported (Khatoon ef al. 2017, 
Shey et al. 2019, Damfo ef al. 2017). 


MATERIALS AND METHODS 

To design multi-epitopic peptide construct, 
immunodominant peptides were preferentially picked 
based on our foregoing laboratory experiments. A series 
of web servers were used to evaluate various interactions 
of the proposed MEV model with potential receptors with 
their physiochemical and immunological parameters. A 
summary of the general workflow is illustrated in the 
schematic diagram (Fig. 1). 


Retrieval of target protein sequences 

In the preliminary experimental set, for designing multi- 
epitope peptides, firstly the amino acid sequences of 
Boophilin (PDB id: 2ODY), (CRPTI_RHIMP) Protease 
inhibitor carrapatin (UniProt id- P81162) and Rms-17 (the 
secretory serpin plays a pivotal role in the tick feeding 
process and pathogen transmission) [NOV41407 (Porto 
Alegre strain, Brazil). 

Supplementarily, to ensure the phylogenetic relatedness 
of the amino acid sequences, target proteins from other 
members of the related tick family were also retrieved. 
To enhance the immunogenicity activation-associated 
secreted protein-1 of Onchocerca 
volvulus (AOA044VDO00), a TLR-2 agonist was chosen 
as an adjuvant. 


Prediction of linear B-cell epitope 

The B-cell epitopes are crucial for vaccine design as 
they stimulate the humoral immune response. Two 
different servers, BepiPred-2.0 (https:// 
services.healthtech.dtu.dk/service.php?BepiPred-2.0) 
and ABCpred (https://webs.iiitd.edu.in/raghava/abcpred/ 
ABC_submission.html) were used to predict B-cell 


Table 1. Presenting statistics of interaction of MAP construct with TLR2 (Cluspro 2.0 server). 


Cluspro 2.0 (Weighted Score) 


Docking with Balanced energy Electrostatic energy Hydrophobic §VdW+Elec 
Receptor (S) score score 
Bovine TLR-2 99 mem -1051 97 members - 1047.9 119 members -1090.1 126 members —275.5 
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epitopes of varying lengths (Jespersen et al. 2017, Saha 
and Raghava 2006). 


Prediction of IFN-y inducing epitope 

Interferon-gamma (IFN-Y) is a pro-inflammatory 
cytokine that plays a potent immunostimulatory role. It is 
accredited as one of the crucial links between adaptive 
and innate immune systems. Apart from increasing the 
level of response to MHC antigens, it is one of the major 
activators for macrophages and Natural Killer (NK) cell 
Prediction of IFN-y epitopes for the designed peptide was 
accomplished using IFN epitope server (http:// 
crdd.osdd.net/raghava/ifnepitope/scan.php) (Dhanda et 
al. 2013). 


Construction of multi-epitope peptide sequence 

From the immunoinformatics prediction, a vaccine 
sequence was constructed using highly conserved 
epitopes. Thereafter, by introducing consecutive GPGPG 
linkers, the different epitopes were combined as a linear 
peptide. Likewise, in an attempt to amplify the 
immunogenicity activation-associated secreted protein-1 
of Onchocerca volvulus, a TLR-2 agonist was integrated 
as an adjuvant (UniProt accession no. AODA044VD0O0) at 
the N-terminus exerting a linker (EAAAK). 


Prediction of immunogenicity and non-allergenic 
property of the construct 

VaxiJen V2.0 (http://www.ddg-pharmfac.net/vaxijen/ 
VaxiJen/VaxiJen.html) tool was used to predict the 
antigenicity of the protein. Antigenicity had been predicted 
by an alignment-free approach based on auto cross- 
covariance (ACC) transformation of protein sequences 
into uniform vectors of principal amino acid properties 
with an input of the final protein sequence as a plain format. 


Table 2. Presenting statistics of interaction of MAP construct 
with TLR2 (HADDOCK 2.4 server). 


Receptor TLR-2 
HADDOCK score -36.14/-9.5 
Cluster size 9 
RMSD from the overall 0.4 +/-0.2 
lowest-energy structure 

Van der Waals energy -35.1+/- 3.6 
Electrostatic energy -329.14/-21.8 
Desolvation energy -7.3 +/- 3.6 
Restraint’s violation energy 721.4+/-72.4 
Buried Surface Area 1577.0 +/- 55.7 
Z-Score -2.2 
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AllerTOP v2.0 (https://www.ddg-pharmfac.net/ 
AllerTOP/method.html) and AllergenFP (https://ddg- 
pharmfac.net/AllergenFP/) were used to predict multi- 
epitope vaccine allergenicity (Dimitrov ef al. 2014a, 
Dimitrov et al. 2014b). 


Evaluation of physicochemical properties 

The physiochemical parameters such as amino acid 
composition, theoretical pI, instability index, in vitro and in 
vivo half-life, aliphatic index, molecular weight, and grand 
average of hydropathicity (GRAVY) were assessed using 
the web server ProtParam (http://web.expasy.org/ 
protparam/) 


Prediction of secondary structure 

The e-services aided by PSIPRED, and Raptor-X 
property were utilized (McGuffin et al. 2000, Kallberg er 
al. 2014). The RaptorXProperty web server (http:// 
raptorx.uchicago.edu/StructurePropertyPred/predict/) 
predicts the secondary structure properties of the protein 
chain through a template-free approach. 


Prediction of tertiary structure 

The input sequence was processed by the RaptorX 
server which enables anticipation of 2D and 3D protein 
structures coupled with DISO and ACC values. 


Refinement of the tertiary structure 

The 3D model of the multi-epitopic peptide was later 
subjected to refinement by GalaxyRefine server (http:// 
galaxy.seoklab.org/cgi-bin/submit.cgi?type=REFINE) 
progressively in two stages. 

Structural changes concerning the initial model are 
expressed in the form of different scores like GDT-HA, 
RMSD, MolProbity. The Clash score, Poor rotamers and 
Rama favored are the three different elements of 
MolProbity score (Heo et al. 2013). 


Validation of tertiary structure 

After selecting 3D structure model, another crucial step 
is the validation of the tertiary structure for the 
identification of any potential error. ProSA-web (https:// 
prosa.services.came.sbg.ac.at/prosa.php) and Molprobity 
server (http://molprobity.biochem.duke.edu) were used for 
tertiary structure validation. The protein structural analysis 
(ProSA web) is a freely accessible tool used for protein 
structure prediction and refinement for theoretical and 
experimental proteins (Wiederstein and Sipp! 2007). 

The molprobity server was used to attain the 
Ramachandran plot for the protein structure (Davis ef 
al. 2007). 
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Prediction of linear and discontinuous B cell a near-native structural computation. Hence, we preferred 
epitope to use two different best-rated servers to validate the 
Linear and conformational epitopes are formed by the results. The basis for using two docking servers was: 
close vicinity of the amino acid residues which were very Haddock 2.4 (https://wenmr.science.uu.nl/haddock2.4/) 
distantly located in the native linear chain, though; owing and ClusPro 2.0. 
to the conformational folding of tertiary protein structure 


they get aligned into the neighborhood. For the prediction Molecular dynamics simulation (MDS) 
of discontinuous B-cell epitopes in the validated protein In an attempt to ascertain the composite atomic motions 
structure, Ellipro (www.iedb.org) web tool was used. of diverse proteins and nucleic acid molecules, a molecular 


dynamics study is done. In this experiment, a molecular 
Molecular docking of the constructs with bovine dynamics simulation study was conducted for the designed 
TLR-2 vaccine construct recruiting the iMODSweb-server (http:/ 
The type and level of molecular interaction of the  /imods.Chaconlab.org/). 
antigenic molecule to that of the immune receptor 
characterize the degree of generated immune responses Codon optimization and in-silico cloning 
and thus identification of active sites of the immune Codon optimization contributes immensely when codon 
receptor is very crucial. In the present study, freely | ¢xpression and usage of the host organism differs 
accessible bio-informatics tools named COPORT and __ significantly from that of the parent organism. It’s evident 
SPPIDER were used for the prediction of active sites in _ that degenerate genetic code provides freedom to multiple 
TLR-2 and construction. For the prediction of the codons of encoding single amino acid; codon adaptation 
aforesaid molecular interactions / probable immunogenic becomes essential for improving the expression rate and 
trigger, structural complex-based docking of ‘designed quality of recombinant protein in E. coli (strain K12). In 
multi epitopic peptide = TLR-=2’ receptors was performed. this experiment we used Java Codon Adaptation Tool 
The docking tools work on different algorithms to derive | (JCat) server (http://www.prodoric.de/JCat) which 
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Fig. 1. Brief summary of general workflow is illustrated in the schematic diagram (Created with biorender.com). 


152 


Exploratory Animal and Medical Research, Vol.12, Issue 2, December, 2022 


at MILF’ Fae Al I VAVYTGYAHCPGOGKLTALERK EI ¥GOQNNEYRSEOLI HOGRLEAWN 


ia co Lo a tr 


aa COoCCCCCCHHCCCECCCCCHHHHHHHHHHHHCCCCeECCCHHHCCCCCEECY 
4s AN GT YMPARGEKENMLEMAWDCKLES SADQAWANQC! FGHSE PRODREGYGENY YT 


aoa a ALARA TAT AOR a RE LPR LY RAP ST RO 


HHH HHCHHCCCEREBRERERERECECCCCERREBERERRCCeCeCcecceccectececcceée 
PTOGOMAWAHTVHI GEGVATGECOEGEATLI WIEH¥Y¥E PeGHMY OEY VoRnoOHREe 


_ im J oe 


ececeececeecacareecaecerecaceCceracteceeceeceeeecaeeedcecececeet 
os KVDR DE YT KR KEL SKE BOLE R KE ARK K PRP eK ECT MPFR YT eRe eR GHKHE BP are 
mo Ee] ra aan rh 


COCHGOCCCCCCOCCOOOCOOOOCO BEE SEC OCOOC CRE EEECOOCOOCOBRBEEES 
GCGGNHDNNTESEBECEGPGPGFI TEGCRGNONATATKEGPSOPGMANMAYT YG 


= me =P oF = 


rift 


iad 
wis 3 
ew COCCCECCECCECCCEEBRERERCCeCCECCCCEBEBeCcececeecece 
& RAGE TORD EGE PGPGTAWVI PTTRS AAT GI PPYW YOR PP OP GP Go 
== = eet aan os 
Legend: 
Si are anti - Conners of botloen 
Pialix can: J CE Caren (A) 
— Cuil Preo: d-saad 
Eel Target Gequence 
HELPITRRALTYAVTCTOCRECKLEALER ER EVGOUHEYRESL TSCELERARE TTA PR CEN RLERAMOCKRL ESA A ORHAR 
ee, ee ee 


f = A a A & al ii oa a "i = tal La J = 


PEGCPCHEPAGOMRECTSCERYTATHEEVEVECLERTAGCTOACKSHHEELPELYEMM PSS ERTHEVACOCYLAFTORAWERTYT 
ee SS i ll —. 


_ = — = 


a 4 cI asa = Pm] ap =] on im ) ie = me mm im 
EI GESVATOCOCSATLIVICHTEPSSRNVEE VIL OAGHFERYORESCTIRRECLERSC LEAR EAAARCPSOPEGECTAFSY SE 
aT ee ee ee ee —— i __ 


a —_——$<$—— = 


DIBEGPOPOACHSH SR ETESERECEODFOPSFETOSCHON ERR TR TEES POPRHARATVYOAR OE TORSO PSPOTAVIFY TE 
Se es = Ee — 


—_=S — a 
Es] = = = ™ Ba] = = = = = = =a Ty om mm 
BAAVGLIPEVTORProrore 
= <= ——— | 
(B] 


Fig. 2. (A) Secondary structure prediction of vaccine constructs using the PSIPRED server having (23.5 % ) alpha-helices, 
(13.23%) beta-strands, and (63.2%) coils. (B) Surface accessibility of the final vaccine prediction using RaptorX Property 
server. The Red color represents the exposed residues (57 % ), and the blue colour represents the buried residues(27 %). 
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provides CAI values (Codon Adaptation Index) and GC 
contents within a range of 0.8-1 and 30-70%, respectively, 
to determine the degree of protein expression (Grote ef 
al. 2005). 

Later, the optimized codon sequence of the designed 
multi-epitope construct was cloned in E. coli plasmid 
vector pET-30a (+). Further, two restriction sites, Ndel 
and HindIII were inserted before (N terminus) and after 
(C-terminus) in the appraised sequence, respectively. 
Finally, the expression in the vector was confirmed using 
SnapGene software (www. Snapgene.com). 


Immune simulation 

C-IMMSIM server (http://kraken.iac.rm.cnr.it/C- 
IMMSIM/index.php?page=1) was used to evaluate the 
immune stimulation and obtain an immune response profile 
by in silico simulation. Vaccine injection without LPS was 
selected and three doses of the vaccine were adjusted at 
45 eight-hour-long steps that equal to a 15 days interval. 


RESULTS AND DISCUSSION 

After meticulously screening a large number of epitopes 
using BepiPred-2.0 web and ABCpred servers, 35 B-cell 
epitopes based on their properties like non-allergenic, 
antigenic, peptide physiobiochemical properties, etc. were 
selected. Finally, three epitopes were selected based on 
their conservancy, exposed nature, and absence of 
transmembrane helices (TMHMM Server v. 2.0). 

IFN-y epitopes were predicted from all the short-listed 
epitopes using the IFN epitope. A total of 5 INF-y inducing 
epitopes were used for the final construct. 

From the above-predicted epitopes, a sequence was 
constructed using linkers. With the addition of GPGPG 
linkers, the different epitopes were combined as linear 
peptides. Three B cell epitopes and 5 INF-y were 
selected to design the construct. After integration of 
adjuvant the final construct was found antigenic (0.706), 
non-allergenic, soluble, and non-toxic. 


Fig. 3. Protein 3D modeling, refinement, and validation. [(A). The 3D model of a multi-epitope vaccine was obtained on the 
Phyre2 server following homology modeling. (B). Refinement: superimposition by the GalaxyRefine server of a refined 3D 
structure (colored) on a ‘crude model’ (gray). (C). Validation: Ramachandran plot analysis showing 90.53% in favored, 6.51% in 
allowed, and 2.96% in disallowed regions of protein residues and (D) ProSA-web, with a Z score of - 4.65]. 
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Fig. 4. Linear B-cell epitopes prediction for the final subunit 
vaccine construct. [Three linear B- cell epitopes have been 
shown in Blue (312-340), Red (266-309) and Magenta color (224- 
246)]. 


The development of multiantigen cross-protective ticks 
vaccine is targeted in some laboratories with variable 
efficacy. Imamura et al (2006) reported that when two 
different serpins (RAS-1 and RAS-2) derived from R. 
appendiculatus were combined, it gave a significant 
protective immunity against ticks. In another vaccination 
trial using a multi-antigen vaccine consisting of four 
salivary antigens (Rm239, Rm76, Rm39, and Rm180) 


Fig. 5. Docked complex of subunit vaccine construct with 
Bovine TLR (PDB ID: 2Z7X: A chain). 


specific antibody production was efficiently stimulated 
against two antigens only, Rm239 and Rm76, conferring 
as the efficacy of 73.2% against R. microplus 
(Maruyama et al. 2017). In an alternative pathway, in the 
ongoing consideration, an effort was made to design a 
multi-epitopic peptide construct comprising 
immunodominant peptides, selected from different salivary 
tick proteins which are preferentially picked based on 
solubility and exposed nature of that protein. Indeed, 
protein abundance is one of the key features for the 
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Fig. 6. The results of molecular dynamics simulation of vaccine construct and TLR-4 docked complex. 
[(A) NMA mobility, (B) deformability, (C) variance (red color indicates individual variances and green color indicates 
cumulative variances), (D) eigenvalues, (E) Bfactor, (F) co-variance map (correlated (red), uncorrelated (white) or 
anti-correlated (blue) motions) and (G) elastic network (darker gray regions indicate stiffer regions) ]. 
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Fig. 7. Schematic diagram of interactions between protein 
chains. [Interacting chains are joined by coloured lines, each 
representing a different type of interaction, as per the key 
above. The area of each circle is proportional to the surface 
area of the corresponding protein chain]. 


selection of antigens, as this type of abundant protein is 
more likely to be processed and presented by MHC 
molecules to the T lymphocytes, and the protein 
abundance correlates with the peptide immunodominance 
(Petitdidier et al. 2019). The RmS-17, the saliva secretory 
serpin is involved in the immune evasion of the tick by 
inhibiting eosinophil, neutrophil-mediated inflammation, 
and platelet aggregation. A probable construct including 
multiple B and INF-yepitopes from the selected tick 
proteins was designed and evaluated various 
physiochemical, and immunological characteristics and 
their interactions with potential receptors using a series 
of web servers. 


156 


(57) Sty! 


Bip (80) 
PReR7I - PepXt - Xhot (155) 
MS: 


= ~ Notl (166) 
“Can wm 


(638) Pwit 


(5858) AsiSI ~ Pyul 


(5432) Seal _ 
(5430) Teper - Xmak ‘sats 
(5215) Nrut 
/ RS) 
Xbal (146)) 
T? promoter 
Sgrat (1574) 
Sent (170) 
(4520) BesSal ~ WARAPE (3838) 
(4356) Peil ~Whal (2255) 
pam aes Bcll* (2269) 
BetOIt (2436) 
(4201) PRPT - Tei sat NreeATIt (2461) 
PapOME (2462) 
| Apal (2466) 
| BesHIT (2666) 
Hal (2761) 
(3362) PpeMt 
(35) Fapl - FapAl Bell (3319) 


Fig. 8. Final vaccine in-silico cloning into the pET30a (+) vector. 
[Te codon sequence of the final vaccine is in red, which is a 
1020 bp gene sequence generated by the JCat server. The 
pET30a (+) expression vector is in black. The codon sequence 
is inserted between HindIII (173) and NdeI (1425), forming a 
clone with a total length of 6285 bp. This image was created by 
SnapGene]. 


The final chimera was improved by employing 
specialized spacer sequences. Previously reported 
GPGPG linker was incorporated between our predicted 
epitopes to produce sequences with minimized junctional 
immunogenicity, therefore allowing the rational design of 
a potent multi-epitope construct. 

It has been reported that tick salivary proteins can 
suppress T cell activities and generate Th2 immune 
responses that may induce immune suppression activity. 
In addition, tick saliva can also suppress interferon (IFN)- 
production by the T lymphocytes (Ramachandra and 
Wikel 1992). So, to combat such immune modulation 
events within a host, we have designed a chimeric vaccine 
incorporating different kinds of epitopes. Immuno- 
informatics analysis of the generated constructs indicated 
that it contains a significant number of high-affinity 
interferon-y (IFN-Y) epitopes and linear B-cell epitopes. 

The molecular weight of the final chimera construct 
was 36.49 kDa and it was predicted to be soluble upon 
expression and was consistent with its simulated 
immunogenicity. The theoretical pI of the predicted protein 
is 9.10 indicating that the construct is acidic. The predicted 
instability index (II) is 23.02 which specifies that the 
protein will be stable upon expression. The computed 
GRAVY and the aliphatic index showed that the protein 
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Fig. 9. In-silico immune simulations in the host in response to MEV administration. [A. Intermittent boosts in B-cell populations 
after secondary and tertiary doses. B. Number of B-cell types in three successive immune responses. C. Helper T-cell populations 
after three shots. D. Cytotoxic T-cell populations after three vaccine administrations. E. Dendritic cells per state. F. Generation of 


cytokines and interleukins at three different stages]. 


contains aliphatic side chains, symbolizing potential 
hydrophobicity. The secondary structure of the final 
construct analysis demonstrated that the protein consisted 
pre- dominantly of coils (63.2%), with 23.5% and 13.23% 
residues being alpha-helices and beta-strands, respectively 
(Fig. 2). In nature, the unfolded protein portion and alpha- 
helical coils have been reported to be one of the important 
forms of “structural antigens” (Corradin et al. 2007). The 
3D structure of the final construct improved strikingly 
after the refinement with the GalaxyRefine server and 
showed desirable results on Ramachandran plot 
predictions. The Ramachandran plot shows that most of 
the residues are in the favored and allowed regions (90.5% 
and 6.5%) with very few residues (2.96 %) in the outlier 
region; a sign of the satisfactory model (Fig.3). 

A recent study has reported the involvement of TLR- 
2 with tick infestation and blood histamine levels in bovine 
(Zhao et al. 2021). Stimulated dendritic cells generate 
lipopolysaccharide, the ligand of TLR-2, and ultimately 
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weaken host immunity predisposing to tick infestation 
(Oliveira et al. 2010). To potentiate immune interaction 
between TLR-2 and the chimeric construct, a data-driven 
protein-protein docking analysis was performed and a 
TLR-2 ligand was used as an adjuvant with the designed 
chimera. Furthermore, molecular docking of the final 
construct with TLR-2 was performed in the ClusPro 2.0 
and Haddock 2.4 servers to compare and validate the 
accuracy, stable protein-protein binding, and molecular 
interactions. It was found that the molecular docking was 
significant with negative energy values (Table-1 and 2, 
Fig. 5) of the top-ordered protein-protein docking complex. 
The lowest energy score of these complexes reflects the 
highest binding affinity between TLR and protein structure 
(Pandey et al. 2018). The molecular dynamic simulation 
results confirmed that the final construct can interact with 
TLR-2 suitably. 

Using the C-ImmSim platform, we evaluated the 
vaccine construct’s capacity to induce a significant 
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immunological response if it was given in a real-world 
circumstance. The findings demonstrated that after the 
primary immune response, the secondary and tertiary 
responses increased steadily. Antibodies (IgM+IgG, 
IgG1+IgG2, and IgG1) concentrations were found to have 
risen significantly (Fig. 9A). Aconsiderable spike in levels 
of antibodies (IgM+ IgG, IgG1+IgG2 and IgG1 antibodies) 
was observed (Fig. 9A). The number of B-cells also 
increased noticeably, and IgG1 biotypes were seen along 
with considerable memory cell development (Fig. 9B). 
Active T-cells were seen to dramatically rise after getting 
a secondary and tertiary injection dose; however, they 
gradually reduced at later stages (Fig. 9C, D). According 
to Fig. 9E, levels of dendritic cells (DCs) and macrophages 
increased significantly following vaccine administration. 
Additionally, increased levels of cytokines (such as (IFN- 
Y) were observed (Fig. 9F). These findings suggest that 
the protein could stimulate a protective immunological 
response. 

Further, the final requirement to authenticate the in- 
silico construct is an expression of the recombinant 
protein in a suitable host. Preferably, Escherichia 
coli (strain K12) expression systems have opted for the 
production of recombinant proteins followed by Codon 
optimization for achieving efficient recombinant protein 
expression within the E. coli expression system. The 
assessment values of the codon adaptability index (0.9731) 
and the GC content (56.14%) are complementing high- 
level protein expression in the E. coli system (Fig. 8). 
However, the experimental validation of the chimeric 
protein is needed to support the present study. A very 
scanty number of such studies have been carried out 
because of the limited availability of validated 
crystallographic protein structures (PDB) belonging to 
different tick species (Aguirre ef al. 2018). To confront 
TTBDs in an efficient and environmentally friendly 
manner, new generation tick vaccine studies should be 
globally endorsed and explored embracing the “One 
Health” approach. 
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